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We derive equations for fluid dynamics from a non-extensive Boltzmann transport equation con- 
sistent with Tsallis' non-extensive entropy formula. We evaluate transport coefficients employing 
the relaxation time approximation and investigate non-extensive effects in leading order dissipative 
phenomena at relativistic energies, like heat conductivity, shear and bulk viscosity. 



I. INTRODUCTION 

The transport properties of matter at extreme high 
temperature and energy density, well proved in high en- 
ergy accelerator experiments at the Relativistic Heavy- 
Ion Collider (RHIC) and at the Large Hadron Collider 
(LHC), raised some fundamental questions related to the 
smallest distance scale in physics. Theories, motivated by 
conformal field theory and higher dimensional dual grav- 
ity, in fact predicted a lower limit for the shear viscosity 
over entropy density ratio being around I/Att in natu- 
ral units Although this limitation can be overcome 
in more sophisticated, nonlinear dual gravity models l2| 
or other approaches the actual calculation [6l-[ll| 

and measurement of the viscosity and related properties 
of elementary matter remained in the focus of research 
interest. 

The more, because the scientific evaluation of experi- 
mental signals, extracted from analyses of particle spec- 
tra and correlations, can be interpreted in terms of ther- 
mal concepts only if our knowledge about the material 
quality of a (strongly) interacting quark-gluon plasma is 
well established [ij]- Besides the interpretation of tem- 
perature [H, [3 , the whole hydrodynamical and statis- 
tical approach relies on our basic assumptions about the 
equation of state (EoS) and the dissipative capacities of 
this new stage of strongly interacting matter (isl . [T6j . 

Modern statistical physics methods to the theoretical 
description of high-energy phenomena include a consider- 
ation of non-extensive thermodynamics [T^-IlEl • Besides 
giving a more general framework than Boltzmann did 
for the entropy formula, also the mesoscopic background 
including Langevin, Fokker-Planck or Boltzmann type 
equations can be extended to involve deviations from 
the classical nineteenth century picture (26l - [30| . Most 
of these generalizations depend on a single parameter, 
denoted by q, which also influences the canonical equi- 
librium state leading from an exponential Gibbs distri- 
bution in energy to a power-law tailed one. These are 
observable in experimental transverse momentum spec- 
tra and multiplicity fluctuations ^3T1 - [33 |. therefore the 
important question arises, how are the kinetic and hy- 
drodynamical approaches affected by the value of this 
parameter. Are there qualitative or only quantitative, 
minor deviations to be expected when q ^ 11 Fits to the 
Relativistic Heavy-Ion Collider (RHIC) and the Large 



Hadron Collider (LHC) spectra make a hadronic matter 
value of qn ~ 1.08.. 1.2 [35-:37j and a quark matter value 
of qg w 1.22 [H] probable. 

These fits point out a deviation from the standard 
Boltzmann-Gibbs (BG) statistics for g = 1. The de- 
gree of non-equilibrium in a system is defined relative 
to a fixed local g-equilibrium. Assuming q ^ 1 means 
that the local equilibrium differs from the standard BG 
equilibrium. On the other hand it may be regarded as 
non-equilibrium compared to the classical q = 1, if one 
considers a time evolution of this parameter. However, 
we do not consider this possibility in the present paper. 

We investigate a general non-equilibrium system which 
- for any given q - dissipates energy and produces entropy. 
This we call g-non-equilibrium. Similarly to classical BG 
systems this g-non-equilibrium is related to the response 
of the system to gradients of different thermodynamical 
intensives. This g-non-equilibrium state does not relax to 
a standard BG equilibrium, but to a local g-equilibrium. 

In this paper we derive the relativistic fluid-dynamical 
equations of motion from a non-extensive relativistic 
Boltzmann transport equation (NEBE) and calculate 
the transport coefficients employing the relaxation time 
approximation for the collision integral. Such calcula- 
tions were done in the non-relativistic limit (soj with 
a non-extensive Boltzmann equation pol l4lj | slightly 
different from ours. The g-generalized first order rel- 
ativistic Navier-Stokes-Fourier equations of relativistic 
dissipative fluid dynamics are also derived within this 
framework. We use the metric of flat space-time 
g^^'^ =diag(l, — 1, — 1, — 1) and natural units throughout 
this work, h = ks = c = 1. 



II. NON-EXTENSIVE BOLTZMANN 
EQUATION AND RELATIVISTIC 
HYDRODYNAMICS 



We follow Lavagno (42| and recall the q-modified non- 
extensive Boltzmann equation (NEBE) compatible with 
Tsallis' suggestion for a generalized non-extensive en- 
tropy formula. Extensions of the classical Boltzmann 
equation, like the one we adopt here from j42j . are put 
forward without any derivation from microscopic quan- 
tum field dynamics, with the purpose to study classical 
statistical effects due to correlations possibly induced by 
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finite system size compared to the characteristic interac- 
tion range jlSj . 

For the sake of simplicity we neglect external field ef- 
fects which simplifies the Vlasov operator in the following 
general Boltzmann-type equation: 



(1) 



where fk = /(x'^ , fc'^) is the single particle phase-space 
distribution function and fk — {fk)'^, with q G (0,2). 
C [/] is the collision integral. The space-time coordi- 
nates are — {t, x) and the particle four-momentum 
= (fc", k) is normalized to the mass m = ^fc^fc^. The 
NEBE is a closed equation assuming that the number of 
binary collisions around space-time coordinates x'^ is pro- 
portional to an expression Hq [fk, fk'], where the explicit 
form of the collision integral involving binary collision 
with incoming momenta p'', p'^^ and outgoing momenta 
fc'', /c'^ is given as, 

C[f] = ^J dK'dPdP' 

X Wkk,-.pp, {Hq [fp, fp'] - Hg [fk,fk']) ■ (2) 

Here the factor 1/2 takes into account the symmetry 
of the initial or final momentum configuration in case 
of identical particles if we interchange the labels i.e., 
o or O {p''',p'')- Further- 

more, dK = (?d^k/ ((27r?i)^fc") is the Lorentz-invariant 
momentum-space volume, with g being the number of 
internal degrees of freedom, e.g. spin degeneracy, while 
Wkki^ppi is the Lorentz-invariant transition rate. The 
transition rate is symmetric with respect to the sequence 
of final states Wkk/-^pp/ = Wkk/^p/p, as well as it is sym- 
metric for time reversed processes Wkk/^ppi — Wppi^kki- 
In this way it fulfils the detailed balance property. The 
collision integral is positive for any / > and it does 
not change sign under time reversal. The g-generalized 
version of the so-called Stosszahlansatz, the assumption 
of molecular chaos, 

Hg [fk, fk'] = exp^ [In, (fk) + In, (fk')] , (3) 

uses the q-deformed exponential and logarithm functions, 

e^Vq{x)^[l + {l-q)x]'/'^'-'\ (4) 

(5) 



"9 
In, {x) 
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For g — > 1 the above expressions are the standard expo- 
nential and natural logarithm functions. Here we stress 
that for g = 1 the detailed balance solution of the above 
equations obey the classical Boltzmann-Gibbs statistics, 
not the usual quantum statistics. 

The entropy four-current has to be properly defined as 



5^ = 



dK fc'' 



A- In, ifk) - fk 



(6) 



demanding positive entropy production according to the 
second law of thermodynamics: 



dK\n,{fk) k^djk 



> 0. 



(7) 



Note that Lavagno defined the entropy four-current as 
S'^ = — J dK kf^fk [In, (fk) — 1], while the current form 
was introduced, correcting Lavagno's definition, by Os- 
ada and Wilk |43l - l46| . We shall return to this point later. 

In order to elucidate the relation to a fluid dynamical 
description, let us deflne now the following momentum- 
space integrals of the NEBE, 



dK k>'djk 



dKC [/] 



(8) 



S^T'^^ = J dK k'^k'^djk = J dKk^C [/] , (9) 

where we introduced the g-modified particle four-current 
and symmetric energy-momentum tensor. 



m= dK k'^fk. 



dK k'^k'^fk. 



(10) 
(11) 



It is straightforward to show that the left hand side of 
Eqs. (HI) and (|9]) vanishes when the particle number, 
the energy and momentum in individual collisions are 
conserved, that is, if pf^ + p'^ = k^ + k''^ holds. This 
property is expressed by specifying the so-called collision 
invariants by using, ijj — a + /3 ^k^^ , such that 



J dKaC[f]^0, 
dKp^k''C[f]=0. 



(12) 
(13) 



It implies the macroscopic conservation laws of particle 
four-current and energy-momentum tensor for any solu- 
tion of the NEBE, 



d^N'' = 0, 



(14) 
(15) 

besides a non- negative entropy production from Eq. ([T]) . 

The entropy production vanishes in local g-equilibrium 
and it leads to the collision invariant ip = In, (/ofe). The 
distribution function which satisfies d^S'^ (/ofc) = is the 
canonical equilibrium one. This result is equivalent with 
the requirement that In, (fok) — cto + P^k^, whence the 
g-equilibrium distribution is given by 



fok = exp (ao - l^ok^'u^) . 



(16) 



Here a = ao and = /3ou'^ are the collisional invariants 
given previously, and it^ is a four-vector normalized to 
one, m'^m^ = 1. These quantities will be identified with 
the inverse temperature /?o = ^/T, the chemical poten- 
tial over temperature, ao = fJ-fT, and the fiuid dynami- 
cal four- velocity of matter u^. The formula in Eq. (H 
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reduces to the well known Jiittner distribution or rela- 
tivistic Maxwell-Boltzmann distribution for g = 1, that 
is fj = exp (ao — Pak^u^). Therefore it is clear that in 
local g-equilibrium there are only five fields which com- 
pletely characterize the system for any given q. 

For the purpose of a quick reference let us introduce 
A^"^ = g^^^ — u'^u^ , used to project an arbitrary four- 
vector into another four- vector orthogonal to u'^. Any 
four- vector, A^^, and in particular the four-momenta of 
particles can be decomposed into two parts using an arbi- 
trary fluid dynamic flow velocity, u'^: k'^ = E^u^ -I- 
where Ek — k^u^ is the Local Rest Frame (LRF) en- 
ergy of the particle and fc^'^^ = A'^'^ki, contains the LRF 
momenta 47[. The LRF or co- moving frame is defined 
with uIj^p = (1,0,0,0). Moreover, the space-time gra- 
dient or four-divergence, 9^ = d/dx^ = u^d/dr + V^, 
should also be decomposed into parts parallel and or- 
thogonal to the flow. Here, the co-moving or proper 
time-derivative, d/dr = is also denoted by an over- 

dot, = dA^ /dr, and = A'^di, denotes the gradi- 
ent. For second-rank tensors the orthogonal and trace- 
less projection is defined as A^'^"^ = A'^''"^Aq^, where 

Now, using the local q-equilibrium distribution func- 
tion we are prepared to introduce the following general- 
ized thermodynamic integrals: 



^ (^2jliy.i J '^^^^^y^'' [A^^'^k^k^y fok, 



(17) 



(2j + 1) 



(18) 
(19) 



where i,j > are natural numbers and (2j-|-l)!! = 
(2j -I- 1)!/ (j!2^) denotes the double factorial. At a 
given fixed temperature {fok^ = (^/ofc/Sao) l^o *^nd 



Ok) 



= {d{fokY /dcto) I A,, whence it follows that 



Jc 



q(i,3) 



fdl, 



/3o 



i.j) 



(20) 



The above integrals for q < 1 run over a finite range 
in energy, Ek < T/(l — g), and are always convergent. 
For q > I the integrands have power-law tails satisfy- 
ing a convergence criterion: (i -I- 2) < l/{q — 1) is re- 
quired for the Iq{i,j) integrals to be finite. Correspond- 
ingly the maximal index i is one higher for the J'g and 
two higher for the JCq integrals. Finally, for g = 1 the 
familiar Jiittner distribution makes all integrals conver- 
gent by its exponential tail for any finite indices. In this 
case = Jq=i{i,j) = ICq=i{i,j) = we re- 

trieve the familiar relativistic thermodynamic integrals, 
see Refs. [48l - l50l |. Note that in the paper by Lavagno 
[4^ [5l| . the <7- modified Bessel function of second kind 



was introduced, Kn{q,z), where z = m/3o, which in our 
notation corresponds to Jq(i,j) in the LRF. 

It is also straightforward to show by partial integration 
that the following recursive relations hold in equilibrium: 



Po 



2i 



^q(i-l,j), (21) 



Po Po 
Furthermore, one can also show that 



^q{h]) ^ (2i + 3)2'g(i+2,i+i): (23) 
(2j + 3) J,(,+2,,+i), (24) 

^q{i+2,j) = "^"^^g(ij) ~ (2j + 3)/Cg(i+2,i+l), (25) 



Jq{i+2,j) = 'm'^Jq{i,j) - (2j + 3) Jg(j+2j + l), (24) 



as well as 



dl, 



— Jq(i,j)O^0 — Jq(i+l,j)Po, 



dl, 



q{i,j) 



(26) 



and similarly 



Jq{^:i) = — ~ "^0 H JTS — Po 



9ao 9/3o 

= K^q{i,j)aa — /Cg(j+ij)/3o. (27) 

Let us recall the definition of the g-modified entropy 
four-current and calculate the equilibrium entropy den- 
sity. So = ^oUfj. with Sq — S*^ ifak) for the g-equilibrium 
distribution function from Eq. (ITCl) . Making use of the 
previously introduced g-generalized thermodynamic inte- 
grals we obtain 



So 



-ao J dKEkfok+f3o J dKElfok + J dKEkfok 



-aoJg(l,0) + PoJq(2M) +lq{l,0)- 



(28) 



Identifying the momentum integrals, or shortly "the mo- 
ments", as the particle density no = 1/5(1,0)7 energy den- 
sity eo = Jq{2,a), and pressure po = -Jq{is)- Here 
we have used Eq. (PT|) to obtain Jq(2,i) = — /^(^^^^(i.o) 
which translates into the familiar ideal gas EoS, po = 
/3(^^Iq(i 0). This is the reason why the definition of the 
entropy four-current has to be given in the specific form 
of Eq. ([5]) . Hence the fundamental thermodynamic equa- 
tion follows directly from Eq. (|28| : 
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-aono -f /3o (eo + Po) ■ 



(29) 



Note that /3oPo 7^ "0, because no = Jg(i,o) 7^ 1-q{i,Q)- 
However, if the number of particles is conserved, the in- 
tegrals of and Jq should lead the same number of 
particles per unit volume, although the normalization of 
the distribution functions differ. 

The definition of entropy leads to the fundamen- 
tal thermodynamic relation and the well known Gibbs- 
Duhem relations, 



dsQ = — aodno -I- l3odeo, 



dpo 



no 
/3o 



dao 



(ep +Po) 
/3o 



d[3o- 



(30) 
(31) 
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The last equation follows from Eqs. (|22l27p . Note that, 
the above thermostatic relations resemble those of classi- 
cal BG thermodynamics. Kinetic theory does not specify 
the zeroth law of thermodynamics, therefore, whether the 
thermodynamic temperature is different from the tem- 
perature in kinetic theory is an open question [52]. How- 
ever, choosing a single parameter for the temperature, 
the equations of g-fluid dynamics become formally iden- 
tical with that of classical fluid dynamics for q = 1. 

Before going further we discuss the fluid dynamical 
equations in g-equilibrium. The particle four-current 
and energy-momentum tensor are calculated from Eqs. 
(jlOlllI) using the local ^-equilibrium distribution func- 
tion from Eq. hence 

= nou^', (32) 
T^'' ^ eou'^u'' - poA^"" . (33) 

These decompositions are formally identical with that of 
a perfect fluid, hence some authors called it a perfect q- 
fluid d^, H^. Moreover, the conservation laws from Eqs. 

[T5|) lead, at least formally, to the well known Euler 
equations of perfect fluid dynamics, 

0, (34) 
•^M^ - n (35) 

making a closed system of equations by supplementing 
with an EoS. These equations imply a vanishing local 
entropy production in g-local equilibrium, 9^5*0 = 
for smooth initial conditions, i.e. without discontinu- 
ities, and EoS's without first order phase transitions. 
Note that in general the EoS is not restricted to be 
an ideal gas EoS. In our case the EoS is given as, 
Po — ^ (eo — '7i^j7g(o,o)) I which in the massless limit re- 
turns the familiar cq — 3po ^^nd hence the speed of sound 
is, Cs = ■\/l/3, irrespective of the q-parameter [42j . 

On this account the general formula can be 
derived rewriting the above equations with the 
help of the thermodynamic integrals and re- 
lations, dao = {deo + ICg^3fi)dPo) /JCg^2fi) and 
d(3o — — (deo — ICq^2,o)dao) //C^(3 o). After some algebra 

we get that the speed of sound squared = (^|f^^ at 

fixed entropy per particle, sq/uq, is given as. 



d^T^ = 0, 



V 



9(1.0) 



3 ^g(2,0) 
-g(3,0)^g(0,0) 



(^«(3,0)^g(0,0) - ^g(2,0)^9(l,0)) 



(36) 



where ho = /C,(3_i)//C,(2,i) = (eo+Po)/fio is the en- 
thalpy per particle and X',(i.o) = K^q{2,o)!<^q{afi) ^ ^Infi)' 



V. 



9(2:0) 



/c, 



9(3,0)^g(l,0) 



N(2,0)- 



fk 7^ /ofe- This will lead to additional terms in the ther- 
modynamic quantities and increase the entropy until the 
system relaxes to local g-equilibrium according to the 
microscopic dynamics described by a NEBE. Once lo- 
cal g-equilibrium is reached these additional dissipative 
quantities vanish. 

Making use of the previously introduced notations, the 
macroscopic fields, such as the particle four-current and 
energy-momentum tensor, can be decomposed in a gen- 
eral frame the following way, 

m = u'^J dKEkh + I dKk'^'^'> fk, (37) 

T^'" = u^'W' J dKElh + ^ A^"^ j dK {A^^k^kfs) fk 

+ 2u^^ f dKEkk'^'''>Uk+ f dRk'^^k'^^fk, (38) 



where the round brackets around the greek indices denote 
symmetrization: = {A^"" + A"'^) /2. Therefore, we 

can uniquely identify the fundamental fluid dynamical 
quantities such as the particle density, energy density 
and isotropic pressure: 



n = u^N^ = / dKEkfk 



e EE u^.u.T^''' = / dKElfk, 



(39) 
(40) 



p = -^A^.T'"- dK (A"^fcafc^) fk. (41) 

Similarly the particle diffusion current, the energy- 
momentum current and the stress tensor are given by 

V^" = A^N" = J dKk^''^ fk, (42) 
W = A>;3T"^ - J dKEkk^^^fk, (43) 
TT^"^ = T^^'^^ = / dRk^^'k^^fk. (44) 



These macroscopic fields, expressed above as momentum 
integrals over the non-equilibrium distribution function, 
fk, are part of the particle four-current and energy mo- 
mentum tensor of a non-equilibrated fluid. Some of these 
quantities are related to their equilibrium thcrmodynam- 
ical pendants through the so called matching conditions. 
It is standard practice to assume that the particle density 
and energy density are unchanged from their equilibrium 
values. 



n = no, e = eg. 



(45) 



III. BEYOND Q-EQUILIBRIUM 

If a system is out-of-equilibrium the distribution func- 
tion is different from the local q-equilibriuni distribution. 



while the isotropic pressure, p, separates into two parts: 
the thermodynamical pressure. 



Po 



^-IJdK {A^^k^kp)fok, 



(46) 
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and the bulk viscous pressure, 



(47) 



such that p = po + n. Here we introduced the devia- 
tion from equilibrium, Sfk — fk — fok- The matching 
conditions require that, / dKEkSfk — J dKE^Sfk — 0, 

while (/fc) = (<5/fc) , W (/fe) = T4^^ ((5/^) and 

TT^" (^fk^ = TTt"" (^Sfk^ . Now, it is clear that only the 

non-equilibrium deviations lead to dissipation, since the 
moments of the equilibrium distribution function lead to 

vanishing dissipative quantities 11 ^/ofc^ = (^fokj ~ 

(/ofe) = (/ofe) = 0. 

For completeness, here we calculate the entropy four- 
current from Eq. ([S]) up to first order in deviations from 
equilibrium <5/fc, by expanding the entropy four-current 
around /qa.-. This expression is given as, 



{Ik) 



differ (ao - iSok^^u^) Sfk + 0{5fl] 



where = —aoN^ + PqTq'^Ui, + (3qPqu^^ and we made 
use of the definitions from Eqs. |32] and |331 Note that 
the higher order deviations from g-equilibrium can only 
be calculated once Sfk is specified. The first order result 
is independent of such details. 

We recall that Eckart [s^ was the first to discuss rel- 
ativistic dissipative fiuids in a proper way. The energy- 
momentum tensor decomposition proposed by him con- 
tained no bulk viscous pressure, 11 = 0, and no parti- 
cle diffusion current, V^^ — 0, hence the fluid dynam- 
ical flow velocity was fixed to the conserved particles, 
= N'^ / ^^/Nt^Nf^. Moreover the energy-momentum 
current in this case was identified with the heat-flow, 
= Later Landau and Lifshitz [53| introduced 

a different decomposition where the flow is tied to the 
flow of energy- momentum = T^'^u^ / ^Tt^'^UaT^pu^ . 
This defines the fiow as the only time-like eigenvector 
of the energy- momentum tensor. Therefore the energy- 
momentum diffusion current vanishes, = 0, in this 
frame. These physically different choices are related by 
the general expression for the heat- fiow = W^ — h[)V^. 

Here we showed that using the NEBE we obtain for- 
mally the same results as in the case of the classical Boltz- 
mann equation, but they are generalized to include the 
non-extensivity of entropy. 



IV. TRANSPORT COEFFICIENTS IN THE 
RELAXATION TIME APPROXIMATION 

Since the equations of dissipative fiuid dynamics are 
not closed due to the fact that fk is unknown, we need 
to find additional equations and relations which close the 
system of equations. Let us assume, following Chapman 



and Enskog [55| . that the non-equilibrium contributions 
from the streaming term are negligible near the equilib- 
rium (5/fc <C /o/c, that is k^^d^Sfk — > 0. In this way 
k^dfj_fk — k^dfj_fok and therefore all non-equilibrium mo- 
ments vanish from the left hand side of the transport 
equation. Hence we obtain the gradients of thermody- 
namic quantities multiplied by transport coefficients. We 
calculate the streaming term in g-equilibrium, 

k^djok = qflr^ {k^d^ao - Ekk^'d^Po - fSok^kf^d^u,) . 

(48) 

From here we follow the pioneering works of Ander- 
son and Witting (56| and employ the relaxation time ap- 
proximation. This simplifies the collision integral of the 
NEBE, 



C [fk] — ^Ek 



{fk 



fok 



(49) 



where tc is the relaxation time. The relaxation time can 
be interpreted as a mean time between collisions. Fur- 
thermore, this approximation fixes the local rest frame 
so that the energy fiow vanishes, as it will be shown be- 
low. Note that these approximations are made in order 
to simplify the collision term, which in the NEBE is even 
more complicated than in the classical case. Neverthe- 
less, these model equations and solutions can be used 
to make very reasonable first estimates for the transport 
coefficients. 

Replacing fk = fok + Sfk into the NEBE and using 
the previous assumptions we get the following transport 
equation. 



Utilizing Eqs. 
Sfk = Tcqff 



k^d^fok = -—Sfk- 
Tc 

and P5|) we arrive at. 



(50) 



2q-l 
Ok 



Ek 3 



[k'^kP - ao + Ekk 



.(51) 



Here we replaced d^Uu — u^u^-|--|^?A^y+(7^,y+w^y where 
9 = V f_iuf^ is the expansion scalar, cr'^'^ — V^^u"^^ shear 
stress tensor and a;'*" — iV^u'^ — Vu'*) /2 vorticity. In 
the last term (3ok^^^ k'^'^ajf^u = 0, due to the fact that 
the vorticity is antisymmetric while k^'^ k^'^ is symmetric. 
Note that we made use of the conservation laws of perfect 
fiuids to calculate the proper time derivatives and hence 
we express them in terms of gradients. Thus, applying 
Eq. (P?)) for different i,j values we obtain 

do = 'T-o25g(2^o) [^^0^^(2,0) - ^g(3,0)] (52) 
/3o = 'n0^q{2.0) [hoK.q{lfl) -~ /C,(2,0)] ^, (53) 

^'^ = Po^ [^^0 'V^ao - V^/3o] . (54) 

Knowing the non-equilibrium distribution function, we 
can calculate now dissipative quantities. The stress ten- 
sor can be obtained in a straightforward manner from 
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Eq. (HH), together with the definitions of the integral 
from Eq. (fTO|) . as being 



2 (rc/3o/C, 



9(3,2)J 



(55) 



whence the traditional Newton-Navier- Stokes relation 
between the stress tensor and shear tensor emerges: 



The shear viscosity coefficient is given by 
V = Tc/3o^q(3,2)- 



(56) 



(57) 



The relativistic generalization of the Newtonian shear 
CT^'', and the Navier-Stokes relation between the stress 
and shear have the same form for the g-modified statis- 
tics as the ones obtained using the classical BG statistics. 

Similarly to the shear viscosity the heat-conductivity 
coefficient can be calculated. One can show that ~ 
(^i?(2,i) ^ ^o^q(3.i)) = in the Anderson- Witting relax- 
ation time-approximation. Using Eq. (|42l) and — 
—hoV^ the result is 



Tc (V^ao) 



..If 



= Tc (V'ao) [/l0^q(l,l) ~ ^g(2,i)] • 
The Fourier-Navier-Stokes law becomes 
9^ = -K (V Vo-'V-ao) 



(58) 



where = -/loT^^ (y^T - Tw^) and the coefficient 

for heat-conductivity or thermal conductivity is given by 



(60) 



Finally the bulk viscous pressure also known as the 
volume or second viscosity, can be calculated from Eq. 
(|47)) where we used Eqs. (|52l53p to express the proper 
time derivatives in Eq. (1ST]) , hence 



n 



-no 



-no 



/3o| j dKflr'Ej:^ {A-^k^kf, 



30 



/10/C2 



'TC-. 



2?20 
5/3o^g(3,2 



^30 ^ ^0^20 



A^20 — hoK. 



10 



20 



3^0 



20 



9(3,1) 



9(2,1) 

(61) 



From this the relativistic version of the classical Stokes 
result follows immediately: 

n = ~C0, (62) 

with the bulk viscosity coefficient, 



no 
3 (/C20 



Jo 



2^20^9(3,2) 



710 

holCia) JC 



3 (/C30 — /l0^2o) ^g(2,l) 
9(3,1)J • (63) 



Once again we note that the above results formally re- 
semble those of the BG statistics, the only explicit q- 
dependence occurs in the definitions of the g-modified 
thermodynamic integrals. All formulas show that the 
transport coefficients are directly proportional to the 
mean free time between collisions and so inversely pro- 
portional to the cross section. 

Having defined the dissipative quantities, the deviation 
from the ^-equilibrium Eq. [51] can be written as, 

= qflr^ (CnH + Cqky^ + C^k^k-^TT^,) , 

where Cn, Cq and are tedious expression inversely 
proportional to the respective transport coefficients, and 
may be calculated from the above formulas. 



V. NON-EXTENSIVE EFFECTS IN THE 
TRANSPORT COEFFICIENTS 

In this section we investigate the previously obtained 
results in the LRF where Ek\LRF — -I- fc^. The 

thermodynamic integrals are calculated in hyperspher- 



ical coordinates, k'^kn = m , k^u^ 



mcoshx. 



(59) fc2 _ m^sinh^Xi A^^'^kfj_k,y = — m^sinh^x and dK ~ 



(27r) 



^47rm^ sinh xd-X- Note that the integrals for g > 1 
are calculated up to infinity while for g < 1 the upper 
limit of integration is Xmax — cosh^^ (^ m(i'-q) ) ■ Further- 
more, the chemical potential is set to zero /i = 0, while 
the spin degeneracy is g = 1. 

In Fig. [1] we plot the natural logarithm of the q-th 
power of the g-equilibrium distribution function from Eq. 
(HH) in the LRF, with m = 140 MeV and T = 140 MeV 
for longitudinal momenta k^ E [0, 10] GeV and trans- 
verse momenta kx = 0. The functions arc plotted for 
q = 0.9 (thick line), q — I (thin line), q — 1.1 (dashed 
line) and q = 1.2 (dotted line). Note that for q = 0.9 the 
function is practically cut-off at fc/, = msinhxmax- 

Next we plot the speed of sound squared as a func- 
tion of temperature and q-parameter from Eq. (j36|) . The 
mass of particles is fixed at m = 140 MeV, while the 
temperature is T € [2, 500] MeV. We see that for any 
value of q the speed of sound is increasing with tempera- 
ture while it saturates for T 3> to to the ultra-relativistic 
limit, Cg = 1/3 independently of q. This is a straightfor- 



ward result from Eq. (p6|. The speed of sound is larger 
for larger q at any given temperature. 
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FIG. 1. (Color online) The fok distribution as function of lon- 
gitudinal momenta for different g-parameters: q — 0.9 (thick 
line), q ~ 1 (thin line), q — 1.1 (dashed line) and q — 1.2 (dot- 
ted line). The mass of particles as well as the temperature 
are m = T = 140 MeV. 
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FIG. 2. (Color online) The as a function of temperature 
for different g-parameters; q = 0.9 (thick line), q = 1 (thin 
line), q — 1.1 (dashed line) and q = 1.2 (dotted line). 



The shear viscosity divided by the relaxation time from 
Eq. (157)) . T]T^^ = /i^o^q(3,2)i is shown in Fig. [3^. On the 
other hand the shear viscosity to entropy ratio divided 

can be in- 



by the relaxation time -r^ = -7f 2<ML 

spected for different values of q in Fig. [St). 

We observe that the viscosity increases with increasing 
temperature for all g-values. This behaviour is more pro- 
nounced for larger g- values, hence there is a monotonic 
increase in the viscosity coefficient with increasing q at 
a given temperature. At the same time it turns out that 
the shear viscosity to entropy ratio at a given temper- 
ature is almost insensitive to different values of q. Fig. 
[3}d demonstrates that for a given relaxation time tc the 
rj/s ratio does not change significantly and at high tem- 
peratures the g-dependence diminishes. This also means 
that the ratio of dissipative to equilibrium quantities such 
as n^'^/{eo+po) ~ {r]/s)/T does not change at high 
temperatures. This is due to the fact that the entropy 
density also increases with increasing q-values by about 
the same factor. In the classical Maxwell-Boltzmann 
massless limit ('7t^^) = PoI(3.2) = 4/ (15eo) and 
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FIG. 3. (Color online) The shear viscosity coefficient divided 
by the relaxation time, jyrj^ (a) and the shear viscosity co- 
efficient over entropy density divided by the relaxation time 
(77/s) Tq^ (b) as function of temperature for different g- values. 



S(q=i) = /3o (^ (9,n ) — -^(2,1)) = 4/3oeo/3 we recover the fa- 
miliar results |57| |. 

Similarly to the viscosity coefficient, we also plot the 
coefficient for the heat-conductivity from Eq. ([SO]) , as 
well as the ratio of heat-conductivity to entropy density 
divided by the relaxation time. Both of these quantities, 
shown in Figs. 2^ and |3)d, are sensitive to the value of 
the parameter q. 

Furthermore, we observe that the heat-conductivity in- 
creases with increasing g-values, and - other than before 
- the ratio -T'p} increases in a noticeable fashion. 

The last figures. Figs. [5^ and Fig. [5Jd show the bulk 
viscosity coefficient divided by the relaxation time, and 
the bulk viscosity coefficient to entropy density ratio di- 
vided by the relaxation time. 

We observe that the value of bulk viscosity coefficient 
is very small. It is about lO^'^ times smaller than the co- 
efficient of shear viscosity C ~ lO^"^??, which corroborates 
the classical result that the bulk viscosity in the non- 
relativistic or ultra-relativistic massless limit vanishes, 
while it is much smaller than the shear viscosity in be- 
tween for massive particles ^57, . i58] . On the other hand 
the bulk viscosity to entropy density ratio does show a 
monotonic increase with increasing q but only for T <^ m. 
For T > m the behaviour is reversed due to the fact 
that the entropy density increases with temperature at a 
slower rate compared to the viscosity. 
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FIG. 4. (Color online) The heat-conductivity coefficient di- 
vided by the relaxation time (a) and the heat-conductivity 
coefficient to entropy density ratio divided by the relaxation 
time as a function of temperature for different g- values. 



VI. CONCLUSIONS AND OUTLOOK 



problems are not solved by introducing the parameter 
q. This manque can be resolved by an appropriate 
choice of the relaxation time in the evolving fluid [60l |. 
Superior but more complicated equations of motion may 
be derived from the NEBE as in the classical theory, 
with methods different from the one presented in this 
work, see for example Refs. j6TI - [65j |. 
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In this paper we have calculated the transport coeffi- 
cients from a g-modified Boltzmann equation. The re- 
sults show that all transport coefficients increase mono- 
tonically with increasing g-values for relativistic ideal 
gases. The speed of sound is also increasing with in- 
creasing q but it saturates quickly to the well known 
ultra-relativistic limit independent of q. 

Furthermore, we derived the g-generalized versions of 
the classical Navier-Stokes-Fourier equations of relativis- 
tic dissipative fluid dynamics from a g-generalized Boltz- 
mann transport equation. These equations relate the dis- 
sipative quantities to the thermodynamical forces linearly 
with positive transport coefficients which dependent on 
parameter q. 

The resulting equations, are nevertheless parabolic, 
hence acausal and unstable [H^. Hyperbolicity related 



FIG. 5. (Color online) The bulk viscosity coefficient divided 
by the relaxation time (a) and the bulk viscosity coefficient 
to entropy density ratio divided by the relaxation time as a 
function of temperature for different g- values. 
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